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Critical exponents in abelian projected SU{2) QCD 

Shun-ichi Kitahara'^,Shinji Ejiri'^,Tsuneo Suzuki^ and Koji Yasuta 
'^Jumonji University, Niiza, Saitama 352, Japan 

''Department of Physics, Kanazawa University, Kanazawa 920-11, Japan 

The critical properties of the abelian Polyakov loop and the Polyakov loop in terms of Dirac string are studied 
in finite temperature abelian projected SU{2) QCD. The critical point and the critical exponents are determined 
from each Polyakov loop in the maximally abelian gauge using the finite-size scaling analyses. Those critical 
points and exponents are in good agreement with those from non-abelian Polyakov loops. 
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Figure 1. Non-abelian, abelian and monopole 
Polyakov loops versus f3 on 16'^ x 4 SU{2) lattice. 



1. Introduction 

Abelian projected QCD is regarded as 
an abelian theory with electric charges and 
monopoles. Dual Meissner effect based on the 
condensation of the monopoles can be consid- 
ered as the mechanism of confinement of quarks. 
This picture is likely to be realized at least in 
the maximally abelian gauge: the value of the 
string tension and the behavior of Polyakov loops 
are reproduced by the abelian link fields and by 
the monopoles (Dirac strings) [|l]-||; the effec- 
tive monopole action is calculated and it indicates 
that QCD is always in the monopole condensed 
phase from the comparison of the energy and the 
entropy 0. 

Figure|l||^ shows the abelian and the monopole 
Polyakov loops defined later appear to be good 



order parameters. However, those curves seem 
to have different slopes. Their absolute values in 
the deconfinement phase are also different. The 
critical behavior of 4-dimensional SU{2) lattice 
gauge theory is shown to be the same as that of 3- 
dimensional Z2 theory. Since the abelian and the 
monopole Polyakov loops seem to be good order 
parameters, it is interesting to evaluate critical 
exponents and critical points exactly from these 
Polyakov loops. What kind of universality class 
will appear from these abelian quantities? 

2. Definitions of Polyakov loops 

After abelian projection is over, we can de- 
fine abelian Polyakov loops written in terms 
of abelian link variables u^{x,t). We separate 
u^{x,t) from gauge-fixed link variables U^{x,t): 
Ui^i{x,t) — A^{x,t)u^{x,t). We define an abelian 
Polyakov loop 

Pabei{xQ) = Re [exp(i^ 6*4(3;, i)J4(a:,i))]. (1) 

Here J^{x,t) — 5i^i,i5x,xa a-nd 9^(x,t) are the an- 
gle variables of u^{x, t) — exp{i9^{x, t)a3), where 
(73 is a Pauli matrix. 

The abelian Polyakov loop can be decomposed 
into two parts: a monopole part and a photon 
part[|j. The abelian field strength ff^i, = d^j^Oi, — 
d^O^ can be separated into two parts: f^i, = f^v+ 
27rn^i^, where n^^ is an integer and /^^ G [— tt, tt). 
Then, rewriting Eq.(|l|), we get 

Pabei{xQ) ^ Re [Pi(xo) • /2(a;o)], 
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Pi(xo) = cx.pi~i^D{x - x',t~t') 
y.d'Mx',t')Ji{x,t)), 

P2{xa) = e^p{~2TTi^D{x~x',t-t') 

xd'^n^A{x' ,t')Ji{x,t)). 

Here D{x,t) is a lattice Coulomb propagator 
which satisfies d'^d^D{x,t) — —^xfi^tfi- The 
monopole Polyakov loop, Pmonoixo) = RePi(a;o) 
is composed of Dirac strings of monopoles. 
Suzuki et al.[^ have indicated that Pabei{xo) 
and Pmono{xo) vanish in the confinement phase, 
whereas Pphoton = RcP2(a;o) is finite at the range 
from (3 —2.1 to 2.5 and did not change drastically 
around the critical point. 

3. Finite-size scaling theory 

We calculated the critical exponent of the non- 
abelian, the abelian and the monopole Polyakov 
loops from a finite-size scaling theory. The sin- 
gular part of the free energy density on x Nt 
lattice has the following form: 

where x = {T — Tc)/Tc. Here the action contains 
the term hNfL [L denotes the magnetization) 
and only the largest irrelevant exponent (— oj) is 
taken into account. By differentiating /s with 
respect to /i at /i = 0, we get 

x{x,N,) = N]l''Q^{xNl'\g^Nr), 
gr{x,Ns) = Qg,XxNl'\g^N^'^), 

where (i), x ^-nd are order parameter, suscep- 
tibility and 4-th cumulant, respectively. Expand- 
ing those equations with respect to x, we have 

i(x = 0,iV,) = 7V7'^/''(co + C3iV7"), (2) 
x(x = 0,iV,) = 7V;/''(co + C3iV7"), (3) 
gr{x = 0,Ns) = gr + c^N~^, (4) 

at a; = 0. The critical point can be defined as 
that point where a fit to the leading A^s-behavior 
has the least minimal ||] . Actually, leading A^^- 
behavior of Eq. (^) and of Eq. (||) is given by 

\nL{x = Q,Ns) = ^{l3/v)\nNs+\nco, (5) 
\nx{x^Q,N,) = (7/z/)lnAr, +lnco, (6) 
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Figure 2. X^ /Nf from g^- fits versus (3 in the 
non-abelian, the abelian and the monopole cases. 
The number of degrees of freedom, TV/ is 2. 

as in ref.|l. From the fits to Eqs.(|), (|) and (|), 
we can find the position of critical point /3c, and 
then obtain the values oi^jvfjjv and g"^ at /3c si- 
multaneously. We also considered the derivatives 
of the observables with respect to x. The lead- 
ing iV^-behavior of each derivatives at the critical 
point is obtained similarly: 

do 

In— (a; = 0,iV,) =plnA^, +lnco. (7) 

Here O is L, x and gr with (p = (1 — 0)1 v, (1 + 
7)/i^ and Ijv). 

4. Results and Discussions 

We performed the numerical calculations on 
Nl X 4 lattices, where TV^ =8,12,16 and 24. The 
standard SU{2) Wilson action was adopted and 
abelian link valuables were defined in maximally 
abelian gauge. We calculated the following ob- 
servables: 

L = ^P(x)/iV3, 

X 

X = N^{{L^)-{m, 
gr = {L^)/{L'f-3, 

where P{x) denotes Psui2)ix), Pabei{x) and 
Pmono{x)- The values of observables at various 
/3 are needed in order to calculate the derivatives 
with respect to x, where x = {f3 — (3c) / Pc- Then, 
we used the density of state method(DSM)|^. 
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1.117(27) 
0.617(16) 
0.311(19) 


1.138(10) 
0.616(12) 
0.299(19) 


1.091(84) 
0.617(57) 
0.326(69) 


1.085(14) 

0.621(6) 

0.326(8) 


1.072(7) 

0.6289(8) 

0.3258(44) 


0.67 
0.35 
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1.977(29) 
3.600(38) 
0.616(25) 
1.218(68) 


2.025(34) 
3.646(44) 
0.617(29) 
1.249(81) 


1.991(88) 
3.608(93) 
0.618(68) 
1.23(19) 


1.944(13) 
3.555(15) 
0.621(8) 
1.207(24) 


1.970(11) 
3.560(11) 
0.6289(8) 
1.239(7) 


1.32 




2.985(47) 


2.995(56) 


3.05(15) 


2.994(21) 


3.006(18) 




V 


1.447(41) 
0.633(13) 


1.438(42) 
0.621(14) 


1.438(41) 
0.600(13) 


1.403(16) 
0.630(11) 


1.41 

0.6289(8) 





Table 1 

The critical exponents calculated from various Polyakov loops. 



First we performed Monte-Carlo simulations at 
/3o=2.2988, and then the expectation values of the 
observables in the vicinity of f3o were obtained us- 
ing DSM. L, X and g,. at /3o were calculated every 
50 sweeps after 2000 thermalization sweeps. The 
number of samples was 100000, except on 24^^ x 4 
lattice (47000 in the case). The errors were de- 
termined according to the Jackknife method di- 
viding the entire sample into 10 blocks (4 blocks 
on 24^ X 4 lattice). 

We estimated the critical point /3c from 
method [s]] . The data of our DSM results are fitted 
to Eq.(^(|) and Eq.(0) at each /3. The number 
of input data is 2 and that of fit parameters is 
2 {u! in Eq.(H) is fixed to 1 in accordance with 
Engels et al.Q). Figure ^ describes the typical 
curves of X^/^f versus /3. 

Averaging the obtained minimal positions of 
X^/Nf, we get 

psu(2) ^ 2.29940(20), 
^abei ^ 2.99962(26), 

prnono ^ 2.29971(23). 

Those critical points are very close to each other. 

Table |^ lists the critical exponents on each crit- 
ical point in the non-abelian, the abelian and the 
monopole case. We get the following results: 

1. The critical exponents in the abelian and 
the monopole case are in agreement with 



non-abelian exponents within the statistical 
error. 

2. Those critical exponents agree with those of 
Z2 rather than those of U{1). 

3. Hyperscaling relations are well satisfied. 

4. Non-abelian exponents obtained are consis- 
tent with those of Engels et al. Q . 

The first and the second results indicate the 
abelian (monopole) dominance in quark confine- 
ment. 
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